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Abstract: We study spin 1/2 isoscalar and isovector, even and odd parity candidates 
for the © + (1540) pentaquark particle using large scale lattice QCD simulations. Previous 
lattice works led to inconclusive results because so far it has not been possible to unam¬ 
biguously identify the known scattering spectrum and tell whether additionally a genuine 
pentaquark state also exists. Here we carry out this analysis using several possible wave 
functions (operators). Linear combinations of those have a good chance of spanning both 
the scattering and pentaquark states. Our operator basis is the largest in the literature, 
and it also includes spatially non-trivial ones with unit orbital angular momentum. The 
cross correlator we compute is 14x14 with 60 non-vanishing elements. We can clearly dis¬ 
tinguish the lowest scattering state(s) in both parity channels up to above the expected 
location of the pentaquark, but we find no trace of the latter. Based on that we conclude 
that there are most probably no pentaquark bound states at our quark masses, correspond¬ 
ing to m,r=400-630 MeV. However, we cannot rule out the existence of a pentaquark state 
at the physical quark mass corresponding to 77^=135 MeV or pentaquarks with a more 
exotic wave function. 
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1. Introduction 

One of the mysteries of hadronic physics has been the failure to observe baryon states with 
quantum numbers that cannot be explained in terms of three quarks. However, for a long 
time this was not considered to be a practical problem due to the presumed large decay 
width of these exotic baryons. The experimental signal of the 0 + (154O) particle [1], [2]- 
[12] changed this situation dramatically. Indeed, the experimental upper bound so far on 
the width of the © + is around 10 MeV. This remarkably narrow width would also explain 
why the 0 + has not been seen before. Since the 0 + was observed to decay into a neutron 
and a K + , its strangeness has to be +1, the third component of its isospin is 0, and its 
minimal quark content is dduus. Prom the lack of a signal in the I 3 = 1 channel the SAPHIR 
collaboration concluded that the 0 + is most probably an isospin singlet state [3]. Its spin 
and parity cannot be pinned down based on currently available experimental data. 

Though the 0 + (154O) is seen experimentally in low energy exclusive processes, there 
are a number of (e + e _ or high-energy proton collision) experiments, where the 0 + (154O) 
is not seen [13] - [20]. The different kinematical and experimental conditions between the 
low energy exclusive experiments (with experimental evidence for 0 + ) and the inclusive 
experiments (usually non-observations) do not allow a direct comparison so that the null 
results do not prove that the positive experiments are wrong [21]. Nevertheless, it is fair to 
say that the experimental situation is not perfectly clear at the moment. Since there are only 
single experimental indications of other exotic pentaquarks (the possible S (1860) state 
reported by the NA49 experiment at CERN [22] and the charmed pentaquark identified by 
the HI experiment at DESY [23]) their existence is even more debated than that of the 
@+(1540). 
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Originally, the experimental search for the 0 + (154O) was largely motivated by the chi¬ 
ral soliton model [24] that predicted for the first time in 1997 a mass of 1530 MeV and a 
width of less than 15 MeV for this exotic S=+l baryon (for an earlier estimate of the mass in 
the soliton approach see [25]). The experimental evidence of the 0 + pentaquark triggered 
a flurry of theoretical speculations about its possible structure, yet unmeasured quantum 
numbers and on the possibility of the existence of other exotic hadrons. A particularly pop¬ 
ular and successful approach is based on different types of quark models [26-28]. Attempts 
have been made to understand the experimental findings by means of baryon-meson bound 
states [30] as well as QCD sum rules [31]. 

These models substantially differ in the properties they predict for the pentaquark state. 
E.g. several models predict positive parity, while other approaches insist on negative parity. 
Clearly, it is of utmost importance to study the 0 + (154O) without any model assumptions, 
based on a first principles non-perturbative approach, i.e. lattice QCD. 

The difficulty of the lattice approach lies in the fact that the 0 + (154O) mass is very 
close to the NK scattering threshold. In lattice QCD one has to use a finite box, implying 
that the continuum of KN scattering states turns into a stack of discrete energy levels with 
the 0 + (154O) embedded somewhere among them. It is then not an easy task to reliably 
distinguish the 0 + (154O) from these nearby scattering states since all the quantum numbers 
coincide. 

There are a few published works on the 0 + (154O) in lattice QCD. Considering the 
difficulties involved, it is not surprising that the results are not in complete agreement. 
Here we collected the main features of these studies, for a more detailed discussion see [39]. 
Except for one, all the lattice studies report a signal in the negative parity channel close 
to the expected location of the pentaquark. Based on the simple fact that the lowest state 
with opposite parity lies much higher, Refs. [32,33] tentatively identify this state with the 
0 + . Others employ finite volume analysis [36] and twisted boundary conditions [35] to 
distinguish between a two-particle and a one-particle state and they conclude that what 
they see is a scattering state. Ref. [35] on the other hand identifies the first excited state 
with negative parity and from its dependence on the volume concludes that it is the 0 + 
resonance. All these works are largely consistent in the lowest masses they find in both 
parity channels, they only differ in their interpretations. 

The only result, which is inconsistent with the rest is that of [34], observing a state 
in the positive parity channel compatible with the 0 + and a much higher state in the 
negative parity channel. We stress that none of the lattice studies so far could identify 
the lowest expected scattering state in both parity channels. This strongly suggests that 
the wave functions, all based on rotationally symmetric quark sources at the origin, do not 
have sufficient overlap with all the low lying states. Another lattice study, Ref. [38], finds 
some evidence that a pentaquark potential based on the diquark-diquark-antiquark picture 
is energetically more favourable than that of the KN picture. All these investigations use 
the quenched approximation, but with different fermion formulations and pentaquark wave 
functions (operators). 

A reliable confirmation of the existence of the 0 + (154O) from lattice studies is achieved 
only if all the states up to above the expected location of the 0 + (154O) have been identified 
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and the © + (1540) can be clearly distinguished from the neighbouring scattering states. It 
is thus clear that a further more comprehensive study is required and this is our aim in the 
present paper. 

Here we use several possible wave functions (operators) that have a good chance of 
spanning both the scattering and pentaquark states. Our operator basis is the largest in 
the literature, the cross correlator we compute is 14x14 with 60 non-vanishing elements. 
In particular we also include displaced, rotationally non-symmetric spatial quark configu¬ 
rations to allow non-zero orbital angular momentum as well as a better separation of the 
scattering states. 

In the positive parity channel the lowest state we can identify is compatible with the 
lowest expected two-particle state and is already significantly above the 0 + . In the negative 
parity channel we can distinguish the two lowest states that both turn out to be compatible 
with the expected scattering states. At the box volumes we use the 0 + is expected to be 
between the two lowest scattering states, but we see no trace of it there. We also carried 
out the analysis for a smaller volume and found that the volume dependence of the energies 
is compatible with all the identified states being two-particle states. 

In conclusion, we identified all the states around the expected location of the 0 + in both 
parity channels and they all turned out to be significantly different from the 0 + . Since our u 
and d quarks were heavier than the physical quarks (corresponding to 777^=400-630 MeV) we 
cannot rule out the possible appearance of a pentaquark state for lighter quarks. Although 
not very likely, it is also possible that a pentaquark state exists with a wave function having 
very small overlap with all our trial wave functions. 

In the present study we chose to work in the quenched approximation again, which 
is known to be quite successful in reproducing mass ratios of stable hadrons [40] - [43]. 
Compared to our previous analysis we improved by three means. In addition to the cross 
correlator technique and the finite volume analysis we increased our statistics by a factor 
of 2-3. 


2. Cross correlators 

In hadron spectroscopy one would like to identify states with given quantum numbers by 
computing the vacuum expectation value of the Euclidean correlation function (O|0(f)0(O)|O) 
of some composite hadronic operator O. The operator O is built out of quark creation and 
annihilation operators. In physical terms the correlator is the amplitude of the “process” of 
creating a complicated hadronic state described by O at time 0 and destroying it at time t. 

After inserting a complete set of eigenstates | i) of the full QCD Hamiltonian the cor¬ 
relation function can be written as 

(O|0(t)O(O)|O) = ^ I <i|0(0)|0) I 2 e -( Ei ~ Eo)t , (2.1) 

i 

where 

0 (t) = e~ Ht 0(0) e m (2.2) 

and Ei are the energy eigenvalues of the Hamiltonian. 
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Note that since we work in Euclidean space-time (the real time coordinate t is replaced 
with —it), the correlators do not oscillate, they rather die out exponentially in imaginary 
time. In particular, after long enough time only the lowest state created by O gives contri¬ 
bution to the correlator. The energy eigenvalue corresponding to that state can be extracted 
from an exponential fit to the large t behaviour of the correlator. 

In principle higher states could also be identified by generalizing the procedure and 
fitting the correlator with a sum of exponentials. In practice, however, that would require 
extremely high precision data, usually not available in lattice simulations. A much more 
realistic solution can be based on the observation that if the operator O happened to have 
negligible overlap with the ground state in the given sector, a single exponential fit would 
yield the first excited state. This, however, is very unlikely to happen by sheer luck, as it 
would require fine tuning. 

It is exactly this fine tuning that can be performed if instead of one operator O one 
considers a linear combination of the form 

n 

n(t) = Y J V i O i {t). (2.3) 

i=l 

The correlator of 1Z can be easily expressed in terms of the n x n correlation matrix 

Cij(t) = 0)) (2.4) 


as 

n 

R(t ) = (nmm = ]r ^c^(t). (2.5) 

i,j =1 


Morningstar and Peardon used this cross correlator to compute glueball masses on the 
lattice [44]. Their procedure was based on the effective mass defined for a general correlator 
as 


m eS 


C{t + At) 
C(t ) 

Let us now consider the effective mass obtained from R(t), 


= ——In 
At 


( 2 . 6 ) 


m(t ) 


'R{t + At)~ 

= -—\n 

Yli,j=0 v i v 'jCij{t + At) 

[ R(t) \ 

At 

Yli,j=0 ViVjC'ij (t) 


(2.7) 


This can be exploited to construct linear combinations that have optimal overlap with the 
ground state or higher excited states. If the correlator contained only n different states, 
the linear combination with the lowest effective mass would yield exactly the ground state. 
In practice this is a good approximation starting already from moderate values of t, since 
higher states die out rapidly. 

A simple computation shows that the stationary points of the effective mass with 
respect to the variables {uj}" =1 are given by the solutions of the generalized eigenvalue 
equation 

n n 

J2Cij(t + At) V j = X'E'CyWvj. (2.8) 

i =1 i= 1 
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Initially we only asked for the lowest effective mass, but this eigenvalue problem can 
have many solutions. It is not hard to interpret them using the following geometric picture. 
Cij(t) and Cij{t + At), being both Hermitian, can be considered to be the components 
of two quadratic forms on the n-dimensional space spanned by the Uj’s. Let us interpret 
Cij(t) as an inner product on this vector space. It can be seen from eq. (2.7) that the 
effective mass does not depend on the normalization of the vector {/ty}, so we can restrict 
it to be of unit length (with respect to the inner product just defined). It is now easy to 
see that the stationary points of the effective mass correspond to the principal axes of the 
second quadratic form, Cij(t + At). In the language of the generalized eigenvalue problem 
this is equivalent to the statement that two quadratic forms can always be simultaneously 
diagonalized in a vector space: there is a basis orthonormal with respect to one quadratic 
form and pointing along the principal axes of the other one. 

Assuming a generic case with no degeneracies the stationary points will have 0,1,2... 
unstable directions and they yield the coefficients of the linear combinations corresponding 
to the ground state and the higher excited states. Of course this statement again is exactly 
true only if there are only n states in the correlator. The importance of corrections coming 
from higher states can be estimated by checking how stable the whole procedure is with 
respect to varying t and At. 

This gives a general method to determine the optimal linear combinations of n operators 
that have the best overlap with the lowest k (k < n) states. The only disadvantage of this 
procedure is that being based on effective masses, it always uses only two points of the 
correlators to extract the optimal linear combinations. On the other hand, once the optimal 
linear combinations have been found the corresponding correlators can be fitted using any 
standard technique. 


3. Details of the simulation 
3.1 Choice of operators 

One of the most important parts of the whole analysis is the proper choice of operators. 
We need a large number of independent operators, which span a large enough subspace 
containing the scattering states and a possible pentaquark state. 

In order to have really independent operators, we used non-trivial wavefunctions for the 
quark fields. The typical operators used in hadron spectroscopy contain quarks at only one 
lattice point with some Gaussian smearing. These operators have automatically zero orbital 
angular momentum and a spin eigenstate can be guaranteed by correctly choosing the Dirac- 
structure of the operator. This, however, gives a very limited set of operators. Moreover, 
some operators, e.g. the one proposed by Jaffe and Wilczek [26] cannot be implemented in 
this way. 

Therefore we decided to use operators, which contain quark fields at different lattice 
sites. In general the five-particle wave function could be any function of the locations of the 
five quarks. However, since the correlation functions are built up from quark propagators, 
we have to restrict ourselves to wave functions, which are products of the individual quark 
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wave functions: 


0 (x 1 ,x 2 ,x 3 ,X4,x 5 ) = qi(xi)q 2 (x 2 )q 3 (x 3 )q 4 (x 4 )q 5 (x 5 ) (3.1) 


Here, for simplicity we omitted the color and Dirac-structure. These are the elementary 
operators for which the correlators can be computed by single Dirac-matrix inversions. A 
general five-quark operator can be written as a linear combination of such elementary 
operators. 

For the individual quark wave-functions we use a simple Gaussian function centered at 
some lattice site: 

qi(xi) = exp 1 °) ) • (3.2) 

It is easy to see that if not all Xi o-s are the same then the operator will not have a 
spherical symmetry and therefore it will create a mixture of angular-momentum eigenstates. 
According to Appendix |j| we can project out angular momentum 1/2 using the projector 

p(Gi)_ 

We had two sets of operators, one with spatially completely symmetric and one with 
antisymmetric operators. Since the cross-correlator of a symmetric and antisymmetric 
operator vanishes we could perform the runs separately for the two sets. It turned out that 
the symmetric/antisymmetric operators had a good overlap with negative/positive parity 
states, respectively. 

Let us simplify our notation further by allowing only quark wave-functions that are 
centered on points of the z axis only. Operators based on such wave-functions have axial 
symmetry and therefore the spin-projection requires a minimal number of extra operators. 
Let 

qi(di,ri,Xi) = exp — ^ ^ ^ , (3.3) 

where z is the unit vector along the z axis. We will usually omit the Xi argument. 

We used the following set of isoscalar operators: 


Oi 

0-2 

o 3 

o 4 

o 5 


e abc [u^(0, 4 ) 0754 ( 0 , 4 )]{u c ( 0 , 4 )s e ( 0 , 4 ) 7 5 d e ( 0 ,4) + (u <-» d)} 
e abc e adx eb9 h [u r (0) 4)C 75 4(0,4)] \u T q (0,4)04(0,4 )}Cs T c (0,4) 


p(Gi) 


abc r T 


U a (0,4)0*754(0,4)]{u c (0,4)s e (N s /2,4)7 5 4(IV s /2,4) + (u ^ d)} 


(3-4) 


p( Gl ) \ e a b c e a d e ebgh[u T ( f, 2 )C^d e (l, 2)][<(-1, 2)Cd h (-l, 2 )}Cs T c (0, 4) 


p{Gi) 


,.abc r n ,T 


4 (0,4)0*754(0,4)] 


x K(0,4) [s e (JV s /4,4) 75 d e (JV fl /4,4) - s e (-N s /4, A)^d e {-N s /4, 4)] + (u ~ d)}] . 


Here C is the charge conjugation operator and the color indices are shown explicitly. 
The first operator is the one used in our previous work [32] with color index contractions 
corresponding to an NK state. 0 2 was introduced in [33]. The third operator is a shifted 
N — K scattering operator with spin projection. The relative displacement of the nucleon 
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and kaon is half of the spatial lattice size N s / 2, so this operator is spatially symmetric. 
The last two operators are the antisymmetric ones. O 4 is based on the proposal [26]. The 
two diquarks are shifted to ±1 from the origin. Finally, O 5 is a shifted N — K operator 
with distance N s / 4. It is first antisymmetrized, then projected to a spin eigenstate. The 
projection of the last three operators requires the computation of 3, 3 and 6 operators, re¬ 
spectively. Therefore we have to compute the correlation matrix of 14 elementary operators 
(except for the elements connecting operators with opposite spatial symmetry). 


3.2 Simulation parameters and results 

We used the standard Wilson gauge action at (5 = 6.0 to generate our configurations. For 
the measurements we used the Wilson fermion action with four different K U) d values for 
the light quarks: 0.1550, 0.1555, 0.1558 and 0.1563. This spans a pion mass range of 400- 
630 MeV. For the strange quark we used a constant n s = 0.1544, which gives the required 
kaon mass in the chiral limit. The lattice size was 24 3 x 60 and for the largest quark mass 
we also performed simulations on a 20 3 x 60 lattice to see the volume dependence of the 
observed states. 

Table [I] shows the statistics we collected 
in the various points. After performing the 
spin and parity projections, we used the di- 
agonalization procedure described in the pre¬ 
vious section to separate the possible states 
in both parity channels. As mentioned ear¬ 
lier the symmetric operators gave a good sig¬ 
nal only in the negative parity channel while 
the antisymmetric operators had reasonable 
overlap only with positive parity states. This 
can be understood since the parity trans¬ 
formation includes a spatial reflection and 
the nucleon-kaon system has a negative in¬ 
ner parity. Therefore we used only the op¬ 
erators 0\ — O3 to extract negative parity 
states and operators O4 — 0 5 for positive parity. 

We varied both t and At required for the diagonalization over a range of 2 — 5 and 
included the systematic uncertainties coming from this variation in the final errorbars. After 
separating the states we had to extract the lowest masses from the individual correlation 
functions. It turned out that for the excited states neither a correlated nor an uncorrelated 
fit with a single exponential (cosh) was satisfactory since in the asymptotic region where a 
one exponential fit could work the data were rather noisy. We used the following technique 
instead. 

If one plots the effective mass log C(t)/C(t + 1) as a function of t , it should show 
a plateau at asymptotically large t values 1 * . It is easy to show that the effective mass 

1 Actually we used a slightly modified ’’effective mass”, namely the solution of the equation cosh(m e g • 

( t — N t /2))/ cosh(m g g ■ (t — N t /2 + 1)) = C(t)/C(t + 1) to get a flat plateau even for t values close to N t /2. 


size 

operators 

Ku,d 

confs 

24 3 x 60 

o u o 2 ,o 3 

0.1550 

242 

24 3 x 60 

04,02,03 

0.1555 

205 

24 3 x 60 

01,02,03 

0.1558 

205 

24 3 x 60 

01,02,03 

0.1563 

205 

20 3 x 60 

01,02,03 

0.1550 

630 

24 3 x 60 

04,05 

0.1550 

250 

24 3 x 60 

04,05 

0.1555 

144 

24 3 x 60 

04,05 

0.1558 

144 

24 3 x 60 

04,05 

0.1563 

144 

20 3 x 60 

04,05 

0.1550 

234 


Table 1 : The collected statistics for the various 
simulation points. 
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approaches its plateau exponentially: 


m ff (i) = m + a ■ exp (—bt) t —* oo, (3.5) 

where m is the lowest mass in the given channel. One can fit the effective masses with 
the above formula and use it to extract the lowest masses. In this way one also uses the 
information stored in the points before the plateau even if the plateau itself is noisy. This 
technique turned out to be very stable and we could start to fit the effective masses at 
* = 2,3. Fig. I illustrates the method for the first two states in the negative parity channel 
for k = 0.1550. 



Figure 1: The effective masses for the first two states in the negative parity channel for k = 0.1550 
and the fitted exponentials. 

In both parity channels we extracted the two lowest masses (which we denote by ?no 
and mi). It is straightforward to define the ratio ct* = rrii/ (uin + mx) which compares the 
possible scattering and pentaquark states to the nucleon-kaon threshold. The experimental 
value of a for the 0 + particle is a e + = 1.07. 

The summary of our results including also the pion, kaon and nucleon masses is given in 
Table |2[ The zero momentum scattering state is just at the threshold. The first scattering 
state with nonzero momentum is expected at 

Ei = yJm? K + 4-7T 2 / (aN s ) 2 + yjm 2 N + Air 2 / (aN s ) 2 . (3.6) 

Its ratio to the threshold is 1.151, 1.166, 1.177, 1.202 for k = 0.1550,0.1555,0.1558 and 
0.1563, respectively for our larger volume. For the smaller volume ( N s = 20) at k = 0.1550 







m n /(m. 


size 

parity 

Ku,d 

am 7 r 

arriK 

am,N 

«o 

aq 

24 3 x 60 

- 

0.1550 

0.296(1) 

0.317(1) 

0.642(6) 

1.01(1) 

1.16(5) 

24 3 x 60 

- 

0.1555 

0.259(1) 

0.301(1) 

0.613(5) 

0.99(1) 

1.16(5) 

24 3 x 60 

- 

0.1558 

0.234(1) 

0.290(1) 

0.592(6) 

0.99(1) 

1.14(8) 

24 3 x 60 

- 

0.1563 

0.185(1) 

0.272(1) 

0.545(7) 

0.98(2) 

1.28(13) 

20 3 x 60 

- 

0.1550 

0.295(1) 

0.316(1) 

0.647(7) 

1.00(1) 

1.24(8) 

24 3 x 60 

+ 

0.1550 

0.295(1) 

0.316(1) 

0.636(6) 

1.16(2) 

1.45(16) 

24 3 x 60 

+ 

0.1555 

0.258(2) 

0.299(2) 

0.615(14) 

1.13(3) 

1.39(15) 

24 3 x 60 

+ 

0.1558 

0.233(2) 

0.288(2) 

0.595(12) 

1.09(5) 

1.39(17) 

24 3 x 60 

+ 

0.1563 

0.184(3) 

0.270(2) 

0.552(10) 

1.14(8) 

1.31(32) 

20 3 x 60 

+ 

0.1550 

0.295(1) 

0.316(1) 

0.647(6) 

1.21(2) 

1.48(12) 


Table 2: The measured pion, kaon and nucleon masses and the ratio of the first two five-quark 
states in both parity channels to the KN threshold. 


this ratio is 1.211. We can see that in all cases the measured mass ratios are consistent with 
the scattering states. The expected and measured volume dependences of the first excited 
state for negative parity and the ground state for positive parity is shown in Fig. |j. 


1 . 
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Figure 2: The volume dependence of the two lowest states in the two parity channels (left panel: 
negative parity; right panel: positive parity). The dashed lines indicate the expected scattering 
states with 0 momentum and the first two non-vanishing momenta. The dotted line shows the 
experimental value of the pentaquark state. 

For the highest quark mass, where we had the largest statistics quark mass we also 
performed the whole analysis for the isovector channel. The extracted masses and their 
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volume dependence turned out to be qualitatively similar to those in the isoscalar channel. 

4. Conclusion 

In this paper we studied spin 1/2 isoscalar and isovector, even and odd parity candidates 
for the 0 + pentaquark using large scale lattice QCD simulations. The analysis needed 
approximately 0.5 Tflopyears of sustained 32 bit operations. 

Before we summarize the results of the different channels one technical remark is in 
order. The @ + pentaquark is expected to be a few % above the NK threshold. Typical 
lattice sizes of a few fermis result in a discrete NK scattering spectrum, with order 10% 
energy difference between the lowest lying states. Thus, we are faced with two problems. 
First of all we have to find the possible pentaquark signal and the nearby scattering states. 
Secondly we have to tell the difference between them. Finding several states could be 
done by multi-parameter fitting 2 or more effectively by spanning a multidimensional wave 
function basis and using a cross correlator technique. The most straightforward way to tell 
the difference between a narrow resonance and a scattering state is to use the fact that the 
former has an energy with quite weak volume dependence, whereas the latter has a definite 
volume dependence, defined by the momenta allowed in a finite system. 

Clearly, any statement on the existence/non-existence or on the quantum numbers of 
the 0 + pentaquark depends crucially on this sort of separation. None of the previous 
lattice investigations on the 0 + pentaquark was able to carry out this analysis. The most 
important goal of the present paper was to do it. 

The individual results in the odd and even parity channels can be summarized as 
follows (based on the statistically most significant, highest quark mass and assuming that 
+ uik) scales with the quark mass). 

1. Odd parity. The two lowest lying states are separated. The lower one is identified 
as the lowest scattering state with appropriate volume dependence (in this case the p=0 
scattering means no volume dependence). This state is 6 a below the 0 + state. The volume 
dependence of the second lowest state is consistent with that of a scattering state with 
non-zero relative momentum. For our larger/smaller volumes this state is 1.8/1.3<r above 
the 0 + state. None of these two states could be interpreted as the 0 + pentaquark 3 . 

2. Even parity. The two lowest lying states are identified. The volumes are chosen 
such that even the lowest lying scattering state is above the expected 0 + pentaquark state. 
The volume dependence of the lowest state suggests that it is a scattering state. For both 
volumes this state is 6cr above the 0 + state. Since the energy of the second lowest state is 
even larger, none of them could be interpreted as the 0 + pentaquark. 

To summarize in both parity channels we identified all the nearby states both below and 
above the expected 0 + state. Having done that no additional resonance state was found. 

2 Note that the energies of the states are determined by the exponential decays of correlation functions; 
extracting several decay rates, which differ just by a few %, from the sum of noisy decays is in practice not 
feasible. 

3 The volume dependence and the larger operator basis of this work suggest, that the odd parity signal of 
our previous analysis [32], which was quite close to the © + pentaquark mass, was most probably a mixture 
of the two lowest lying scattering states. 
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This is an indication that in our wave function basis no 0 + pentaquark exists (though it 
might appear in an even larger, more exotic basis, with smaller dynamical quark masses or 
approaching the continuum limit). 


A. Parity projection 


In this appendix we summarize how the energy of the lowest state can be extracted sepa¬ 
rately in the two parity channels. Although we consider spin 1/2 baryon correlators here, 
our discussion can be generalized to other states. 

The basic objects one can compute on the lattice are Euclidean correlators of the form 
(O a (x)Og(0)) corresponding to the amplitude of the process of creating a state at time 
zero with the operator Og( 0), evolving it to a later time .To and annihilating it with O a (x). 

There are two complications when one wants to extract the lowest state in a given 
parity channel. Firstly, simple baryonic operators usually couple to both parities, therefore 
one has to project out parity by hand. Secondly, the box has a finite time extent T 
with (anti)periodic boundary condition. Therefore, a single source at time zero is in fact 
mathematically equivalent to the sum of an infinite number of identical sources located 
at t = 0,±T, ±2T.... Due to the exponential fall-off of correlations, only the two sources 
closest to the sink, i.e. at t = 0, T give appreciable contributions to the infinite sum. If we 
assume, as we shall always, that 0 < xo < T then 

OO 

^2 (O a (x,x 0 )Op(0,nT)) « (O a (x,x 0 )Op( 0,0)) + £bc( O a (x, x o )6p(0, T) ), (A.l) 

n =—oo 


where £bc is +1 for periodic and —1 for anti-periodic boundary condition in the time 
direction. The first term on the r.h.s. represents particles propagating from time 0 to xq 
while the second term represents antiparticles propagating from time xq to T. Thus even 
after projecting to a given parity channel the correlator has contributions not only from 
particles of that parity, but also from the antiparticles of particles of the opposite parity. 
Therefore, an additional “projection” is needed to get rid of the latter. 

Before starting to describe in detail how the two projections can be carried out let us 
discuss the form of the first term on the r.h.s. of eq. ( A.l ). We shall assume that O a is a 
spin 1/2 baryon operator and a is its Dirac index. Due to the transformation properties of 
O a , the most general form the correlator can have is 


{O a (x)Op( 0))= [f(x 2 )x^ + g{x 2 )l ] a/3 , (A.2) 

where / and g are scalar functions of the length of the four-vector x and 1 is the 4x4 unit 
matrix. After projection to the zero momentum sector this becomes 


where 


C a p{x 0 ) = d 3 x (O a (x)Op(0)) = [A(x 0 )'y 0 +B(x 0 )l] afS , (A.3) 


A(x 0 ) = J d 3 x f{x 2 )x 0 , B(xq) = J d' ! x g(x 2 ). (A.4) 
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The important point here is that upon integration over 3-space the terms in the correlator 
proportional to the spacelike 7 matrices vanished due to their antisymmetry. We also note 
that from eq. ( |A.4| ) A(xo) and B(x 0 ) are easily seen to be anti-symmetric and symmetric 
respectively. 

Parity projection. Parity projection of an arbitrary operator O can be performed by 

V ± 0 = i [O ± POP" 1 ) , (A.5) 

where P is the parity transformation and V±0 couples to parity eigenstates of ± parity. 
In particular, on spin 1/2 fermionic operators parity acts as 

PO a (x 0 ,x.)P ~ 1 =ri('yo) a vO u (x 0 ,-X-), (A. 6 ) 


where 7 = ±1 is the internal parity of O depending on the parity convention we choose for 
the elementary fields and on how O is constructed from those. The analogous formula can 
be easily obtained for O a . 

When constructing correlators it is enough to project to a given parity either at the sink 
or at the source. This by itself ensures that only states of the given parity are propagating 
in the correlator. Inserting the projection into the correlator of eq. (A.3) at the sink one 
obtains e.g. for the positive parity correlator 

\ (1 + 770 ) [A(x 0 ) 7 o + -B(x 0 )l] = | [A(x 0 ) + t]B(xo)] 1 + ^ [A(x 0 ) + 7 -B(x 0 )] 70 - (A.7) 

The negative parity channel can be constructed analogously by replacing 7 with —7 every¬ 
where. 

Notice that all the matrix elements of the parity projected correlator have the same 
functional dependence, 1/2(A ± rjB), on xo- The exponential fit to this function will yield 
the lowest state in the given parity channel. In practice the simplest way to obtain the 
parity projected correlator is to compute two suitable elements of the 4x4 correlation 
matrix C aj g(x 0 ) that yield A(x 0 ) and B{x 0 ) respectively. So far we pretended that the box 
is infinite in the time direction and neglected the second term on the r.h.s. of eq. (A.l). 

‘Particle” projection: Now the full correlator, including the term that “comes back” 
through the time boundary, the only object that we can actually compute in a finite box, 
has the form 


C(x 0 ) + £b c C(x 0 - T) = [A(x 0 ) - £bc A(T - x 0 )] 7 o + [B(x 0 ) + £b c B{T - x 0 )] 1, (A. 8 ) 

where we arranged the arguments of A and B to be non-negative using their (anti¬ 
symmetry. 

Were we to use our prescription above for parity projection, i.e. the “70 component” of 
the correlator ± 7 X its “1 component”, we would end up with the parity projected correlator 

^ [A(x 0 ) + rjB(x 0 )\ + [—A(T — x 0 ) + B{T — x 0 )], (A.9) 

which, due to an extra minus sign, does not have the simple functional form f(xo)+f(T— xq) 
that could be fitted with a cosh. This extra sign, however, can be easily canceled if we 
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compute the A and the B components (the one proportional to 70 and 1 , respectively in 
eq. ( |A.3| )) with opposite boundary conditions. In that case the parity projected correlator 
has the form 

^[A(x 0 ) + r]B(x 0 )] + i [A(T - x 0 ) + rjB(T - x 0 )]. (A.10) 

If 1/2[A(xo)+7]B(xq)\ is a sum of exponentials corresponding to the energies of the states in 
the given channel, (|A.10| ) is a sum of cosh’s with the same exponents. This is the functional 
form we have to use for fitting when extracting masses. 

B. Projection to a spin eigenstate 

In this appendix we outline how a specific spin eigenstate can be projected out from a given 
lattice hadron operator. After summarizing the relevant group theoretical principles we 
discuss how our spin 1/2 pentaquark operators were constructed. 

When discussing spin on a hypercubic lattice the first problem is that due to the absence 
of full 50(3) rotational symmetry it is not straightforward to assign spin to a lattice energy 
eigenstate. States on the lattice can be classified into irreducible representations of the cubic 
group O or its double cover 2 0, not 50(3) and SU( 2 ) as in the continuum. 

With the exception of the lowest four representations, when restricted to 2 0, irre¬ 
ducible representations of 50(2) do not remain irreducible. The spin 0, 1/2, 1 and 3/2 
50(2) representations are the exceptions, these restricted to 2 0 are equivalent to the ir¬ 
reducible representations A\,Gi,T\ and H, respectively. Also any state belonging to an 
irreducible representation of 2 0 has components belonging to several different spin repre¬ 
sentations of 50(2). For instance a state in G\ has components in spin 1/2, 7/2, 9/2... 
50(2) representations and H has components of spin 3/2, 5/2, 7/2.... 

This means e.g. that if on the lattice we find the lowest energy state in the G\ repre¬ 
sentation of 2 0 , we can identify that with a spin 1/2 state in the continuum, provided all 
the higher spin states contributing to G\, i.e. s = 7/2,9/2... can be assumed to have much 
higher energy. In this sense, for practical purposes, the lowest few representations of 50(2) 
and 2 0 can be identified as follows: 

0<->i4i, 1/2 « Gi, 1 > Tj, 3/2 <-> H. (B.l) 

The task we have at hand is thus to construct states belonging to specific representations 
of the cubic group 2 0. This can be most easily done using the technique of projection 
operators that we summarize here for completeness. The simple form of the method of 
projectors we present here can be used only when each irreducible representation occurs 
in the decomposition at most once. Therefore it is essential to know ahead of time the 
irreducible representations occurring in a tensor product and their multiplicities. This can 
be most easily found using group characters. See e.g. [46] for explicit formulae and character 
tables of O and 2 0 . 

Let G be a finite group, D 2 - ^ ( g ) be the matrix elements of its irreducible representation r 
of dimension d r . Let the transformations T(g ) form an arbitrary (not necessarily irreducible) 
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unitary representation of G. We would like to project a specific irreducible representation 
r out of the carrier space of the T(g)’s. Let us define the transformations 

t 1 = ]§f E T <»). 


(B-2) 


g&G 


where |G| is the number of elements G has and * denotes complex conjugation. 

It is straightforward to show that if \i/j) is any vector belonging to the carrier space of 
T(g )’s then for a fixed j the d r vectors 

i = 1, ...d r (B.3) 


l^)=4 r) l 


either transform as basis vectors of the irreducible representation r or they are all zero. For 


the proof see any standard text on group representations, e.g. Ref. [45]. Equations (B.2) 
and ( |B.3| ) can be exploited to project out different representations of 2 0 from a given state 
on the lattice and its rotated copies. 

In particular, we would like to construct pentaquark states belonging to G± that cor¬ 
responds to spin 1/2. Although more complicated cases can also be considered, here we 
restrict ourselves to the one where the spin indices of all the quarks but one have been 
contracted to be scalars and the total spin of the pentaquark arises by combining the spin 
1/2 (Gi) of the remaining quark with the orbital angular momentum of all the constituents. 
Therefore we have to project G i out of G\ ® s , where s is a representation of the cubic 
group O (not 2 Ol), corresponding to the orbital part. 

In practice s depends on the spatial arrangement of quark sources and this can be 
exploited to make things as simple as possible. Eq. ( |B.2j ) implies that, in general, projection 
to a specific spin involves as many terms as the number of elements of the group 2 0, i.e. 48. 


The situation, however, is much better if the projection formula (R3) is applied to a state, 
with an orbital part having some degree of symmetry under cubic rotations. The simplest 
case is when the five quark sources all have complete rotational symmetry, i.e. the orbital 
part is trivially s = A\. Then all the rotated copies of the quark sources are identical, the 
sum in eq. ( |B.2| ) can be explicitly computed and the projection reduces to projection to 
spin up or spin down. The decomposition here is A\®G\ = G\. All the operators used in 
lattice pentaquark spectroscopy so far fall into this category. 

To explore the possibility of non-zero orbital angular momentum we have to consider 
less symmetric quark sources. Another possibility is to put the antiquark at the origin 
with a rotationally symmetric wave function, displace the two pairs of ( ud) quarks along a 
coordinate axis (say z) keeping the arrangement cylindrically symmetric with respect to the 
z axis. Inspired by the Jaffe-Wilczek diquark-diquark-antiquark picture [26], in anticipation 
of orbital angular momentum 1, we construct this state to be antisymmetric with respect 
to the interchange of the two displaced quark pairs. Let us call such a state | ± z). It is 
easy to see that the rotated copies of this state span a three dimensional space carrying the 
representation T\ of O. A possible set of basis states is given by (ud) pairs displaced along 
the three coordinate axes; | ± x), | 4= y), | ± z). This arrangement corresponds to projecting 
out the spin 1/2 (G*i) component from the decomposition 


Ti ® Gi = Gi © H. 


(B.4) 


14 




(s~1 \ 

Let us choose \ip) = | |) ® | ± z) and compute P,, 1 \il>). The transformations T{g) appear¬ 


ing in eq. (B + ) are direct products of G\ transformations acting on the quark spin and 
transformations acting on the orbital part. Each term in the sum and as a consequence 
the whole sum itself can be decomposed into three terms proportional to | ± x), | ± y) and 
| ± z). The G\ matrices can be easily obtained by restricting the defining representation 
of 51/(2) to 2 0 and with the factors D^' ,> *(g) they can be summed independently for the 
three terms resulting in 


df' 1 [lt)®l±2>] = L ojl T >®l=M+ (°o) 1 T >®I ±y> 


+ J" IT)®I±,> 


= 11) © | ± x) + i\ |) © | ± y) + | T) © | ± z). (B.5) 

In a similar fashion we obtain the other (spin down) basis element of the G\ projection; 

p 2i l] [IT)®I ±2 )] = IT)®|±®> - *|T)®|±J/> - |i)®|±«). (B.6) 

Note that, up to some numerical factors coming from the normalization of spherical har¬ 
monics, these expressions are identical to the spin 1/2 part of the SU( 2) Clebsch-Gordan 
decomposition 1 © 1/2 = 1/2 © 3/2. We could also construct spin 1/2 from similar, but 
symmetric orbital states for the quarks displaced to x,y,z = —d. This would correspond 
to A\®G\ = G\ or 0 © 1/2 = 1/2 for SU{2). 

Building these states requires seven quark sources; an antiquark at the origin and six 
quark sources, two along each coordinate axis (we use the same mass and source for the u and 


d quarks). Eq. (ELI) shows that keeping the same spatial arrangement the representation 


H corresponding to spin 3/2 could also be projected out. However, we have not explored 

(G W 

this possibility here. For that we would have had to replace the matrix elements D\- 1 ( g) 
with those of H. 


m eq. 

Besides the diquark-diquark-antiquark wave function we also wanted to study triquark 
quark-antiquark states. The simplest non-trivial way to do that is to displace the quark- 
antiquark pair along a coordinate axis, say +z. Let us call the orbital part of this state 
1 + 2 ;). Its rotated copies span the six dimensional space with a possible basis formed by 
|+x), |— x), | +y), | —y), \+z), | —z). This space, however, can be split into an antisymmetric 
part spanned by combinations of the form | + x) — | — x), ... and a symmetric one spanned 
by | + x) + | — x), etc. The representation of O on the antisymmetric part is Ti in exactly 
the same way as in the diquark-diquark case, resulting again in the spin projected state 

p 1+ [lt)®(l + z}-|-2»] = 

U) ® (I + x) - I - x)) + i\ |) © (| + y) - | - y)) + | T) ® (| + z) - I - z)). (B.7) 

The three dimensional symmetric part of the orbital space is reducible to A\®E. Thus 
we can also produce spin 1/2 trivially from the symmetric part by A\ © Gi = G\, 


P‘S t) [|T)®(l+z> + |-z»] = 

I T) ® [| + x) + I - x) + I + y) + I - y) + I + z) + I - z)\ . 


(B.8) 
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